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This paper discusses free carrier generation by pulsed laser fields as a mechanism to switch the 
optical properties of semiconductor photonic crystals and bulk semiconductors on an ultrafast time 
scale. Requirements are set for the switching magnitude, the time-scale, the induced absorption as 
well as the spatial homogeneity, in particular for silicon at A = 1550 nm. Using a nonlinear absorption 
model, we calculate carrier depth profiles and define a homogeneity length 4om- Homogeneity length 
contours are visualized in a plane spanned by the linear and two-photon absorption coefficients. Such 
a generalized homogeneity plot allows us to find optimum switching conditions at pump frequencies 
near !//c=5000 cm -1 (A = 2000 nm). We discuss the effect of scattering in photonic crystals on the 
homogeneity. We experimentally demonstrate a 10% refractive index switch in bulk silicon within 
230 fs with a lateral homogeneity of more than 30 /im. Our results are relevant for switching of 
modulators in absence of photonic crystals. 

PACS numbers: 42.70.Qs, 42.65.Pc, 42.79.-c 



INTRODUCTION 



There is a fast growing interest in photonic crys- 
tals; composite materials whose refractive index 
varies periodically on length scales that match 
the wavelength of light. [l| The optical proper- 
ties of photonic materials are determined by the 
spatially varying refractive index, analogous to 
the periodic potential for an electron in a crys- 
tal. Large spatial variations of the refractive index 
cause a strong interaction between light and the 
composite structure. Bragg diffraction causes the 
photonic dispersion to organize into bands, much 
like the energy levels of electrons in semiconduc- 
tors. A major goal of the field is the realization 
of three-dimensional (3D) structures that possess 
a photonic band gap.0, 0] At frequencies inside 
the band gap, the optical density of states van- 
ishes. This should completely inhibit spontaneous 
emission of sources inside the photonic crystal. 
Indeed, strong modifications of the spontaneous 
emission lifetime of quantum dots have recently 
been demonstrated in photonic crystals. Q| In the 
presence of weak controlled disorder^ Anderson lo- 
calization of light is also predicted. Q In this case, 
a photon may be trapped at a point defect which 
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serves as a cavity with a high quality factor. 

In the examples above, the photonic crystals 
themselves do not change in time. Switching ex- 
periments, in which the properties of photonic 
crystals are modified on an ultrafast timescale al- 
low many interesting new opportunities. Switching 
3D photonic crystals is particularly interesting, as 
it provides the dynamic control over the density 
of states inside the crystal as well as a change in 
Bragg reflections. |6( Ultrafast control of the den- 
sity of states should allow the switching of spon- 
taneous emission of light sources inside a crystal, 
and capturing and releasing light in cavities inside 
the crystal. In absence of photonic crystals, ul- 
trafast switching of bulk semiconductors finds ap- 
plications in high speed optical modulators |jj and 
waveguides, Last but not least, ultrafast control 
of photonic crystals is important for controlling the 
propagation of light, such as in switched macrop- 
orous silicon, or 2D crystal slabs. [T(i| 

In optical switching experiments, four important 
requirements have to be met.[(| First of all, the 
magnitude of the induced change in the real part 
of the refractive index n' must be large enough to 
obtain the desired effect. A relative change in n' of 
5% is required to induce a major change in the den- 
sity of states. Such large changes can be induced 
by free carrier generation but are not achiev- 
able with Kerr switching. The second requirement, 
which is important for applications, is the mini- 
mum time scale At over which the switch occurs. 
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In experiments where light pulses are trapped in- 
side photonic crystals, a switching time scale on 
the order of a few hundred femtoseconds is neces- 
sary. The third requirement is that the absorption 
of probe light, gauged by n", should be small in 
the switched sample. Excited carriers in a semi- 
conductor bring about inevitable absorption (re- 
lated by a Kramers-Kronig relation to the change 
in n'), that should remain within limits by a lim- 
ited carrier density. The fourth requirement con- 
cerns the spatial homogeneity of the change An' in 
a sample. Homogeneity is particularly important 
in switching the density of states in 3D photonic 
crystals. A large gradient in An' in the crystal re- 
sults in a highly chirped switched sample, which 
can no longer be considered a photonic crystal. 
All of these requirements also pertain to other ap- 
plications of switched semiconductors, such as in 
waveguiding,Q albeit much relaxed. Therefore, we 
expect the results from the present study also to be 
relevant for applications outside photonic crystals. 

Switching experiments in photonic crystals were 
pioneered by Leonard et a/..[£| Optical free carrier 
generation was used to change the refractive index 
of the backbone of a 2D silicon photonic crystal. 
The carrier density generated in their experiment 
was sufficiently high to induce a shift of a Bragg 
stopgap on ultrafast time scales. Nevertheless, the 
experiments revealed a serious limitation: the ab- 
sorption of the pump beam limits the volume of 
switched material and leads to a spatial inhomo- 
geneity in the degree of switching. Leonard et al. 
deduced that only a layer of three unit cells near 
the sample surface was switched. 9] Similar inho- 
mogeneity is probably also playing a role in recent 
studies of silicon infiltrated opaline crystals, where 
a disappearance of the Bragg peak was observed in 
the absorption range. |12j Therefore, the use of two- 
photon absorption was proposed as a way to in- 
crease the penetration depth of pump light into the 
sample and improve the switching homogeneity. 

In this paper, we investigate the spatial homo- 
geneity of optically generated free carrier plasmas 
in semiconductors. We discuss a non-linear ab- 
sorption model that takes into account both linear 
and two-photon absorption processes. From this 
model, we derive optimum pumping conditions and 
we define a homogeneity length scale to obtain ho- 
mogeneous switching conditions. In particular we 
trace constant homogeneity lengths in a general 
parameter-diagram that pertains to any semicon- 
ductor. We discuss the role of disorder-induced 



diffusion of the pump beam in photonic crystals. 
Finally we experimentally demonstrate that suffi- 
cient ultrafast refractive index changes can be ob- 
tained with sufficient lateral homogeneity. 

We concentrate on silicon, not only because of its 
wide technological use, but also because it allows 
greater homogeneity than GaAs that was discussed 
in Ref. 0- Examples are given for light in the tele- 
com band at ^t c i c /c = 6450 cm -1 (A = 1550 nm), 
but can easily be generalized to other frequencies. 
We set a minimum volume of 5 unit cells cubed 
in which the change An' must remain within 10% 
of its maximum value. This homogeneity require- 
ment holds for both the lateral directions x and y, 
as well as for the z-direction defined in figure ^ 
Even in such a small crystal volume, the DOS al- 
ready shows a significant decrease for frequencies 
that lie in the band gap.[L| The typical unit cell 
size of silicon inverse opal photonic crystals with 
a bandgap near f te i e is a = 1.2 /xm, therefore the 
homogeneously switched area must extend at least 
6 jjim. in all three dimensions. 



II. FREE CARRIER GENERATION 

In optical free carrier generation a pump pulse 
is absorbed by a semiconductor sample, creating 
a free carrier plasma with electron-hole density 
N c h . The generated carrier plasma changes the di- 
electric function e{uo) of the sample by an amount 
Ae c h(w). In the case of silicon, the Drude model 
gives a excellent description of the resulting c(uj) 
for densities below 10 22 cm~ 3 (see Ref. ITTI) . 



e(w) = e B {io) + Ae eh (uj) = e B - ( — ) — i— , 

\ lo / 1 + i— — 

(1) 

where e B is the bulk dielectric constant, u) 
{=2~kv) the frequency of the probe light, m e the 
electron mass, m* pt the optical effective mass 
of the carriers, td the Drude damping time.[T4| 

and uj p = \J {N e he 2 ) / {to m* opt m e ) the plasma fre- 
quency. Under the condition that 1/(o;td) <C 1, 
which is valid for silicon at Wteic, we can derive the 
following simplified expression for the refractive in- 
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dex from Eq. QJ: 



absorption length l ca of the excited carrier plasma 
is equal to 



I'+in" = V^R = Je B -(^) 2 +i - r 



e B -(^) 2 
(2) 



1 



2k n"(w) 



(4) 



For silicon with carrier densities below 
10 20 cm -3 , the refractive index is linear with 
the carrier density within 0.2%: 



£ s - - ; — ; ?N eh . (3) 

Thus, the induced change in the refractive in- 
dex n' (via e(u>)) is completely determined by the 
optically induced carrier density. 

An example of a carrier induced change of re- 
fractive index in bulk silicon is given in Fig. [2] In 
this experiment, a powerful ultrashort pump pulse 
was focussed to a spot with radius r pump = 70 /im, 
resulting in a peak intensity at the sample inter- 
face of Io = 115 GWcm -2 . The reflectivity of a 
weaker probe pulse with a smaller spot radius of 
r P robe = 20 /im was measured in the center of the 
pumped spot at different time delays with respect 
to the pump pulse. The scan in Fig. [3] shows that 
the reflectivity of the sample changes from 32% 
to 28%. The 10% - 90% rise time is 230 fs, con- 
firming an ultrafast change in n'. From Fresnel's 
formula we find the refractive index change to be 
more than 10%, corresponding to a large generated 
carrier density of 1.6 x 10 20 cm~ 3 . This is about 
twice the carrier density that is needed to obtain a 
change in n! of 5% for probe light at Lo tc \ a . We have 
shown that optical carrier generation can be used 
to change n' by a large amount on a sub picosecond 
timescale, meeting the first two of the four main 
density of states switching requirements. 

For applications such as optical modulators and 
waveguides, much smaller changes in n' are already 
sufficient, typically An' = 10~ 4 , see Ref.0 This 
corresponds to a carrier density of 10 16 cm -3 . For 
such low carrier densities, the required pump pulse 
energy at f pnmp /c — 12500 cm -1 is on the order 
of several nJ, allowing the use of diode lasers with 
repetition rates exceeding 1 GHz as a pump source. 
Thus it seems that carrier-induced optical switch- 
ing may have much broader applications beyond 
photonic crystals. 

Next, we will discuss the induced absorption in 
switched semiconductors, which is the third den- 
sity of states switching requirement. The carrier 



where fco = u)n'(uS)/c. After inserting the ex- 
pressions for n'(u)) and n"{oj) from Eq. we get: 



T D C. 



(5) 



We immediately see that the absorption length 
is inversely proportional to the plasma frequency 
ujp squared and thus inversely proportional to the 
carrier density. In silicon, the carrier absorption 
length of probe light at wtoic is ^ca = 22 /jm for 
a carrier density of 10 20 cm -3 . In our analysis 
of 3D photonic crystals we assume the crystals to 
be inverse opals with a typical filling silicon fill- 
ing fraction of $ = 25%. The carrier absorption 
length in such crystals will then be approximately 
four times larger than the bulk carrier absorption 
length, or £ ca — 88 /im. This result shows that 
in our analysis, where the carrier density is less 
than 10 20 cm~ 3 , the carrier induced absorption for 
Wteie remains very small. We have now shown that 
the first three requirements for successful switch- 
ing experiments can be met for silicon at telecom 
frequencies. The remainder of this paper will dis- 
cuss the fourth requirement: the spatial switching 
homogeneity. 



III. OPTICAL PROPERTIES OF SILICON 

To analyze the switching homogeneity in the z 
direction, we briefly consider how pump light is 
absorbed in semiconductors. At low pump inten- 
sities I, the absorption of light in semiconductors 
scales with the intensity: oil, where the absorp- 
tion coefficient a tends to zero for photon energies 
Tiuj below the electronic bandgap energy E gap , see 
Fig.[3]for silicon. At high pumping intensities, non 
linear two-photon absorption starts to play an im- 
portant role. For two-photon absorption, the ab- 
sorption is proportional to f3I 2 where (3 is the two- 
photon absorption coefficient. This coefficient is 
expected to vanish for Tilq > E gap /2. Fig.|21displays 
the frequency dependence of the linear and two- 
photon absorption coefficients of silicon. The data 
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was obtained from our measurements and from 
Refs.Hi 13 13 El: With our lasersystem,[ijj 
z-scan measurements |2(J were performed to ob- 
tain the two-photon absorption coefficient (3 at 
two additional wavelengths. For bulk silicon at 
fpump/c = 5000 cm" 1 , (3 was measured to be 
0.20 ± 0.05 cmGW" 1 . At u pump /c = 6250 cm" 1 , 
we obtained (3 = 0.80 ± 0.1 cmGW" 1 . The latter 
value is in excellent correspondence with the value 
(3 = 0.88 ± 0.13 cmGW" 1 at u pump /c= 6494 cm" 1 
from Ref. Il6l Our measurements confirm that 
f3 tends to zero for photon energies approaching 
Eg a p/2. Both the linear and the nonlinear absorp- 
tion coefficient can thus be controlled by varying 
the pump frequency. 



IV. HOMOGENEITY OF SWITCHED 
SEMICONDUCTORS 

A. Homogeneity in the z-direction 

We now present a model that calculates the car- 
rier density depth profile N e h(z) caused by absorp- 
tion of pump light. The absorbed intensity is de- 
scribed by the nonlinear differential equation 



dl(z) 
dz 



-[al(z)+f31 2 (z)}, 



(6) 



which we have solved by implicit integration. 
The resulting expression describing the intensity 
depth profile is 



I(z) = 



z e 



l + (f3I /a)(l-e—)- 



(7) 



where Io is the intensity at the interface. The 
resulting carrier density profile N c h(z) is related to 
the intensity depth profile I(z) as 



N ch (z) = ^ z)Tpump 



a 



(8) 



where r pump is the pump pulse duration. The 
factor 1/2 for two-photon absorption indicates that 
two photons must be absorbed to generate one 
electron-hole pair. Substitution of Eq. (0 into 
Eq. (JHJ) allows us to calculate the carrier density 
profile for any given combination of a, /? and Iq. 



We now investigate the relation between the 
absorption coefficients and pumping homogeneity. 
To quantify the homogeneity of a switched sample, 
we first define the homogeneity length £hom within 
which the carrier density remains within 10% of its 
surface value: 



0.1 x 



1 dN ch (z) 
_Neh(z) dz 



(9) 



z=0 



Because the homogeneity length is directly re- 
lated to the maximum gradient in the carrier depth 
profile, it is a helpful parameter in quantifying 
the homogeneity of switched semiconductors. For 
applications where switching homogeneity is im- 
portant, the homogeneity length should be much 
larger than the size of the switched sample. To il- 
lustrate the homogeneity length, three carrier den- 
sity depth profiles are shown in Fig. 0] for which 
(3 is kept constant at 2 cmGW -1 , and a is var- 
ied from zero (dotted curve) to 400 cm" 1 (dashed 
curve) to 800 cm" 1 (solid curve). For each case, 
Io is chosen such that the carrier density reaches 
0.9 x 10 20 cm" 3 at the sample interface, corre- 
sponding to a 5% change in n' at v te i e in silicon. On 
the right y-axis the corresponding real part of the 
refractive index n' for silicon at f to i c is shown. The 
pump frequency was assumed to be 5000 cm" 1 . 
A closer look at the three depth profiles in Fig. 0] 
shows that for a = cm" 1 , the homogeneity length 
is 1.0 fim. As a increases to 400 cm -1 , the homo- 
geneity length increases to 1.2 fim. As a increase 
further to 800 cm" 1 , the homogeneity length de- 
creases again to 1.0 /im. The surprising occurrence 
of a maximum in the homogeneity length can be 
explained with the aid of Eq. J8J). For small a 
(a <C Io/?), the absorption is dominated by two- 
photon absorption. If a increases, the pump inten- 
sity Io needed to obtain the surface carrier density 
decreases, reducing the slope of the intensity pro- 
file at the interface, determined by the exponent 
— (a + Io[3). This leads to an increase in 4om- If ol 
increases further to the regime where the absorp- 
tion is dominated by linear absorption (a 3> Io/?) , 
any further increase in a will result in a decrease of 
the homogeneity length. In the region between the 
two extremes, the homogeneity length apparently 
attains a maximum value. This means that sim- 
ply choosing two-photon (or perhaps even higher- 
photon) absorption over linear absorption is not 
always sufficient to ensure an optimal homogene- 
ity. 
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We make a homogeneity plot to obtain further 
insight in the influence of a and /3 on the homo- 
geneity. First we choose a fixed electron density at 
the interface N oh (0) = 0.9 x 10 20 cm" 3 . The cor- 
responding homogeneity length contours are then 
deduced from our absorption model, and visual- 
ized in a plane spanned by linear and two-photon 
absorption coefficients, see Fig. [5] The absorp- 
tion coefficients for silicon at various frequencies 
taken from Fig. |3| are also plotted in the plane in 
Fig. [5] The three depth profiles shown in Fig. 0] 
correspond to positions a, b, and c in Fig. [SJ To 
obtain a certain minimum homogeneity, the ab- 
sorption coefficients must remain below the curve 
corresponding to the particular minimum homo- 
geneity. This graph thus allows us to directly ob- 
tain the homogeneity length that can be obtained 
for a semiconductor at a certain frequency. The ho- 
mogeneity plot demonstrates how pumping homo- 
geneity can be maximized by choosing the appro- 
priate pump frequency. Generally, smaller absorp- 
tion coefficients lead to an increased homogeneity. 
This increase comes at the price of a higher neces- 
sary pump intensity Io- From Fig. [S] we conclude 
that the most homogeneous switch for silicon can 
be achieved for i/ pump /c = 5000 cm -1 . For this 
pump frequency, the homogeneity length is 2.9 fim. 
For comparison: if v pump /c is equal to 1 2500 cm , 
the homogeneity length is only 0.6 (im. The neces- 
sary pump intensity remains below the maximum 
available pump energy of our laser system I max - 
Our generalized homogeneity plot is valid for all 
materials and allows us to find optimum switch- 
ing conditions. The choice for pumping frequencies 
which are low in the two-photon absorption regime 
drastically increases the pumping homogeneity. 

We now make a plot that shows surface carrier 
density contours (N o h(0)) corresponding to a par- 
ticular fixed homogeneity length. From such a plot 
the maximum homogeneous change in refractive 
index can be derived. We choose a fixed homo- 
geneity length of f.5 /im, since inside a typical 
photonic crystal with 25% filling fraction, the ho- 
mogeneity length will be four times larger (6 /zm) 
thus fulfilling our homogeneity requirement. Fig.HJ 
displays contours for two different carrier densities 
oh = 1 x f0 20 cm" 3 and 2 x fO 20 cur 3 . 

In Fig. the upper curve is the homo- 
geneity contour, and the lower one of each pair 
of curves indicates the minimum absorption co- 
efficients for which carrier density can be gen- 
erated given the maximum available intensity 



Imax = 1 TWcm~ 2 .pjj For feasible switching ex- 
periments, the absorption coefficients must be in 
the area to the right of the intersection of the two 
curves. As the carrier density increases, the con- 
stant ^hom hue moves towards lower absorption co- 
efficients, while the I max line moves towards higher 
values. As an example: for a carrier density of 
2 x f0 20 cm~ 3 , homogeneous switching can only 
be achieved for absorption coefficients within the 
shaded area in Fig.El It is seen that this area does 
not overlap with the trajectory of silicon parame- 
ters. Therefore, this carrier density is not achiev- 
able given 4om = 1.5 /im. With decreasing carrier 
density, the range between the curves will over- 
lap the silicon parameter trajectory at some point. 
Such an intersection determines the upper limit to 
the carrier density (given 4om = 1.5 /im) as well as 
the pump frequency that pertains to the relevant 
(a,/3)-point. For silicon, this intersection is calcu- 
lated to occur at N h = 1.9 x fO 20 cm~ 3 . From 
Eq. 10 we obtain the corresponding maximum ho- 
mogeneous change in refractive index at vtc\o to be 
11%. 

The condition of a maximum intensity I max can 
be relaxed by choosing a smaller pump spot ra- 
dius r pump . The assumed pump radius of 75 /im 
provides a lateral homogeneity which is large com- 
pared to the homogeneity requirement, therefore 
we could choose a smaller r pump while maintain- 
ing sufficient lateral homogeneity. A higher pump 
intensity would allow homogeneous switching ex- 
periments at even lower absorption coefficients, al- 
lowing larger changes of the refractive index. How- 
ever, we must keep in mind that the carrier ab- 
sorption length is inversely proportional to the in- 
duced carrier density, see Eq. JSJ. For a refractive 
index change of 11%, we predict the carrier ab- 
sorption length inside a photonic crystal to drop 
to £ ca = 38 /im, which may be sufficient to meet 
our third requirement for switching experiments in 
photonic crystals. At higher carrier densities how- 
ever, the carrier absorption length will become too 
small to meet this requirement. 

A similar analysis can be done for other semicon- 
ductor materials. The carrier density needed for a 
5% change in n' will depend on material proper- 
ties such as the optical effective mass of the carriers 
and the Drude damping time td (see Eq. (QJ). We 
briefly discuss switching of GaAs at a frequency 
v/c — 9430 cm^ 1 proposed in Ref. 0. For GaAs 
the optical effective mass m* opt = 0.06 is 2.5 times 
smaller than for silicon. Therefore, the carrier 
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density in GaAs is 2.5 times lower than what is re- 
quired for the same refractive index change in sili- 
con. We find that the area increases in which the 
homogeneity requirement can be met. However, 
the two-photon absorption coefficient for GaAs 
at v/c = 9430 cm -1 is 26 cmGW" 1 , far above 
the homogeneous switching area.[2l] This imme- 
diately shows, that GaAs at this pump frequency 
cannot be used in homogeneous switching experi- 
ments. From this brief analysis we learn that the 
pumping frequencies in GaAs should be reduced, 
to where the two-photon absorption coefficient is 
much lower value. 



B. Lateral homogeneity 



We now consider the homogeneity in the lat- 
eral (x,y) directions. We have measured the rel- 
ative change in probe reflectivity from a bulk sili- 
con sample pumped in the two-photon absorption 
regime at i/ pllmp /c = 6250 cm -1 . For this fre- 
quency, (i = 0.8 ± 0.1 cmGW" 1 and a = 0. The 
probe delay is fixed at 10 ps to avoid transient 
effects. The pump focus position is shifted with 
respect to the probe focus by adjusting the pump 
mirror with a micrometer drive. The pump energy 
Epump is 10.8 /iJ on a focus with r pump = 80 /im. 
The probe frequency ^ pr obc/c is 7692 cm -1 , with a 
focus size of r pro be = 25 /im. The dashed line in the 
graph is drawn to indicate the width of the pump 
intensity distribution. The measured relative re- 
flectivity shows a minimum at the center of the 
pump distribution, and decreases away from the 
center. The reflectivity data was fitted by a Gaus- 
sian curve (solid curve). The radius of this mea- 
sured reflectivity minimum is 34 ± 5 /im, which is 
considerably smaller than the r pump . To obtain the 
lateral homogeneity length, we have determined 
the maximum lateral distance by the pump focus 
center for which AR/R remains within 10% of its 
maximum value (dotted lines). This distance turns 
out to be 16 /jm in both directions, corresponding 
to a homogeneity length of 32 /im. This illustrates 
that for a sample pumped by a pump beam with 
r pump = 80 /im, the lateral homogeneity length is 
much better than the 6 /im which we required for 
switching of the density of states. 



V. SWITCHING HOMOGENEITY IN 
REAL PHOTONIC CRYSTALS 

The above analysis was done for bulk semicon- 
ductor samples. To obtain the homogeneity length 
in photonic crystals, the bulk absorption length 
is divided by the semiconductor filling fraction 
$. In case of a 5% switch of the refractive in- 
dex in a silicon photonic crystal with <E> = 25% 
at t'pump/c = 5000 cm -1 , the homogeneity length 
thus increases from 2.9 /im to 11.6 /im, which is 
twice the homogeneity requirement of 6 /im. For 
a higher pump frequency of 12500 cm -1 , we find 
a homogeneity length of 2.3 /im, which is too low. 
This illustrates that the homogeneity required for 
switching of the density of states in silicon photonic 
crystals can only be performed at two-photon ab- 
sorption frequencies. 

In the analysis so far, the extinction of pump 
light due to random scattering inside the photonic 
crystals was neglected. We now discuss how to 
incorporate inevitable scattering in photonic sam- 
ples. Scattering is quantified by the mean free path 
£ m f p : the characteristic length over which a co- 
herent beam becomes diffuse. The homogeneity 
length of light inside a photonic crystal is related 
to £ mip and 4bs as 



4om = 0.1 x 



1 

^abs 

weak 



1 



L-mfp 



(10) 



In the limit of weak scattering, where 
4if P > 4bs, EqUUI reduces to 4 om = 0.1x4 bs . In 
the limit of strongly scatter, where £ m f p <C 4bsj the 
homogeneity length becomes i + horn = 0.1x£ m f p . 
As opposed to the adverse effect of scattering on 
the homogeneity in the z-direction, scattering will 
generally be favorable for lateral homogeneity, as 
pump light which is removed from the coherent 
pump beam is be scattered laterally. 

Recently, our group has developed a quantita- 
tive model of the mean free path inside photonic 
crystals. [23 One of the main results of the anal- 
ysis is that the mean free path for pump light 
decreases with pump frequency squared a;~ u 2 mp . 
The absolute value of the mean free path depends 
not only on frequency, but also on many proper- 
ties of the crystal such as the unit cell size, the 
amount of disorder and the refractive index con- 
trast. As an example, we apply the model to a 
silicon inverse opal photonic crystal with a lat- 
tice parameter of a = 1240 nm, corresponding 
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to a bandgap frequency near Wtele- We assume 
combined size poly dispersity and lattice displace- 
ments of 2% (which is beyond the current state 
of the art). From the model of Ref. I22I we ob- 
tain a mean free path of 5.3 unit cells for pump 
light at j/ pllmp /c = 12500 cm -1 . The correspond- 
ing homogeneity length, dominated by £ m f p , would 
thus be 0.5 unit cells, well below our homogene- 
ity requirement of 5 unit cells. For a 2.5 times 
lower pump frequency of 5000 cm -1 , the model 
predicts a (2.5) 2 times larger mean free path of 
^mfp = 33.3 unit cells inside the crystal. To- 
gether with the earlier obtained absorption length 
of 4tbs = 116 fim (= 93 unit cells), we obtain a ho- 
mogeneity length of [(l/93+l/33)- 1 ]/ 1 = 2.4 unit 
cells. This result points out that even at pump 
frequencies near the two-photon absorption edge 
of silicon, the homogeneity requirement of 5 unit 
cells can not be met in silicon inverse opals with a 
lattice parameter of a = 1240 nm. 

The scattering model predicts that decreasing 
the unit cell size will result in higher homogeneity, 
as this will reduce the relative pump frequency. 
Therefore, we consider silicon inverse opals with 
a reduced lattice parameter of a = 900 nm, cor- 
responding to a bandgap near 8900 cm" 1 , just 
below the absorption edge of silicon. The scat- 
tering model predicts a larger mean free path for 
pump light with a frequency of 5000 cm -1 of 
^mfp = 63 unit cells. Using the absorption length of 
116 /im (= 116 unit cells), we obtain homogeneity 
length of: 4om = [(1/116 + l/eS^/lO = 4.1 unit 
cells, close to our homogeneity requirement. 

Major improvements in switching homogene- 
ity can be made by studying diamond-like 
structures. [23l |24| as the lattice parameter can be 
as low as 600 nm for a bandgap near w t cic in such 
crystals. In such structures, the scattering model 
remains unchanged, apart from a constant pref- 
actor that depends on the shape of the unit cell. 
We predict that reduced lattice parameters will 
decrease the relative pump frequency sufficiently 
to allow homogeneous photonic density of states 
switching experiments. 

VI. CONCLUSIONS 

We have discussed four important requirements 
for free-carrier induced optical changes in semicon- 
ductors, with emphasis on density of states switch- 
ing in photonic crystals: the amplitude of change 



in refractive index, the timescalc on which the 
switch takes place, the induced absorption, and 
the homogeneity of the induced change. We have 
demonstrated that the first two requirements can 
be met: a 10% change in n' in bulk silicon observed 
within 230 fs. We have also experimentally demon- 
strated sufficient lateral homogeneity in a bulk sil- 
icon. Furthermore, we have measured two-photon 
absorption coefficients bulk silicon. The induced 
absorption was deduced to be low for carrier den- 
sities below 10 20 cm~ 3 . 

We have discussed a non linear absorption model 
to describe the spatial homogeneity of optically 
generated electron-hole plasmas in semiconduc- 
tors. We have introduced a homogeneity plot, 
which directly relates linear and two-photon ab- 
sorption coefficients to the maximum homogeneity 
that can be achieved for any semiconductor. From 
such a plot, we conclude that for density of states 
switching in silicon photonic crystals, the optimum 
carrier density is about 10 20 cm -3 . To obtain 
the required homogeneity, the absorption coeffi- 
cients must be minimized, by a judicious choice of 
pump frequency. Due to peak intensity limitations 
the lowest pump frequency that can be chosen is 
around 5000 cm -1 . 

We have discussed the effect of scattering in pho- 
tonic crystals on the pumping homogeneity The 
homogeneity of switched photonic crystals turns 
out to be limited by scattering. We conclude 
that the homogeneity condition can barely be met 
in silicon inverse opals. We predict that in dia- 
mond structures the relative pump frequency will 
be small enough to allow homogeneous density of 
states switching experiments. 

Finally, we have briefly discussed the ramifica- 
tions for waveguides and modulators and we con- 
clude that free-carrier switching is also useful for 
applications outside photonic crystal. 
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VIII. FIGURES 



FIG. 1: Schematic image of a sample: The z-axis is 
defined to be pointing into the sample, the x- and in- 
directions are the lateral directions. In the figure, the 
front face of the sample displays a SEM image of a dry 
etched Si 2D photonic crystal with a slightly rhomboid 
symmetry (lattice angle 85 deg). The lattice parame- 
ter a is equal to 750 nm. SEM courtesy of L.Woldering. 
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FIG. 2: Time resolved reflectivity measurement on 
bulk Si, pumped at !/ p „ mp /c = 12500 cm -1 , pulse en- 
ergy E pump = 2 pij, pulse duration r pump = 120 fs, 
r P ump = 70 fim and peak intensity 115 GW/cm 2 
(upper panel). The reflectivity of a probe with 
fprobc/c = 7692 cm -1 , r probo = 20 ^m and 
T P robo=120 fs decreases from 32% to 28%, corre- 
sponding to a calculated carrier density N c h — 1.6 x 
10 cm at the surface of the sample (see right-hand 
scale). The time difference between 10% and 90% of 
the total change is 230 fs. The lower panel shows the 
intensity autocorrelation function (ACF) of the pump 
pulses. The full width half maximum (FWHM) of 
200 fs corresponds to a pulse width of 140 fs FWHM. 
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FIG. 3: Absorption coefficients a and fi of Si versus 
pump frequency. The dashed vertical lines correspond 
to photon energies of E gap /2 and E gap . The open 
squares indicate linear absorption coefficients taken 
from Ref. llRlf right-hand scale). The solid triangles in- 
dicate the two-photon absorption coefficients that we 
have determined by z-scan measurements (left-hand 
scale). The left-filled triangles indicate the values from 
Ref. Uri the right-filled triangles are data from Ref. UtI 
and the open triangle data from Ref. llH 
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FIG. 4: Calculated carrier density depth profile N e h(z) 
for constant (3 = 2 cmGW -1 and three different 
linear absorption coefficients: a = (solid curve), 
a = 400 cm -1 (dashed curve) and a = 800 cm -1 (dot- 
ted curve). With N e h at the interface kept constant 
at 0.9 x 10 20 cm -3 , the necessary pump intensity 
It) is calculated, assuming that r pump = 120 fs and 
r pump = 75 /im. The corresponding carrier density 
depth profile was obtained with Eq. 0. The homo- 
geneity length £hom is shown for a = 0. The right- 
hand scale shows the resulting refractive index n' for 
v/c = 6450 cm -1 light in Si. 
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FIG. 5: Homogeneity plot. Contours of constant 
4om are plotted in the (a, (5) plane, denning regions 
in which homogeneous switching can be achieved. 
The generated carrier density is kept constant at 
N eh (0) = 0.9 x 10 20 cm -3 , enough for a 5% change in 
n' . The solid curve corresponds to ^hom = 2.9 /xm, the 
dashed curve to 4om = 1.0 /im and the dotted curve to 
4om = 0.6 fim for bulk silicon. The closed squares con- 
nected by the dotted arrow are linear and two-photon 
coefficients for Si obtained from Fig. [5] The open cir- 
cles a, b and c correspond to the depth profiles N e h(z) 
plotted in Fig. 0] Pumping parameters: r pump = 120 
fs, r pump = 75 /j,ra and v pU mp/c = 5000 cm - . 
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FIG. 6: Carrier density contours for a particular bulk 
homogeneity chosen as £hom = 1-5 /im. The con- 
tours correspond to two different values for N c h(0): 
the dashed lines to N e h(0) = 1 X 10 20 cm" 3 , and the 
solid lines to N ch (0) = 2 x 10 20 cm -3 . The upper one 
of each pair of curves corresponds to £hom = 1.5 fim. 
The lower curve of each pair indicates the minimum 
absorption coefficients for which the given N o h(0) can 
be obtained without exceeding the maximum inten- 
sity I m ax = 1 TWcm -2 . The closed squares con- 
nected by the dotted curve line are the a and 
coefficients for Si obtained from Fig. El Pumping 
parameters: r pump = 120 fs, r pump = 75 /jm and 
Vmp/c = 5000 cm -1 . 
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FIG. 7: Relative probe reflectivity measurement 
on a bulk Si sample for different lateral positions 
of the probe focus with respect to the probe focus. 
The probe delay is fixed at 10 ps. Pumping pa- 
rameters: Tpump = 120 fs, fpump/c = 6250 cm" 1 , 
Epump = 10.8 /iJ, r pump = 80 /ttm, r pump = 120 fs, 
^probo/c = 7692 cm -1 , r prob c = 25 /im. The squares 
are the measured datapoints, fitted with a Gaussian 
curve width of 68 ± 10 /im (solid curve). The hori- 
zontal dotted line indicates the level where the change 
in reflectivity has decreased by 10%. From the verti- 
cal dotted lines, which indicate the intersections of the 
90% line with the Gaussian fit of the data, we obtain 
a lateral homogeneity length of 32 pirn, illustrating the 
excellent lateral homogeneity. The dashed curve indi- 
cates the measured width of the pump focus. 



